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We obtain the QCD quark condensate from consideration of unquenched quark dynamics in 
Dyson-Schwinger gluon vacuum. We consider the non-local extension of the condensate and deter- 
mine the quark virtuality. We also obtain the condensate-driven contribution of the non-perturbative 
QCD to Euler-Heisenberg Lagrangian of QED in external electromagnetic fields. 
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I. OVERVIEW 

We describe a method relying on Dyson-Schwinger 
equations (DSE) for obtaining (non-local) quark con- 
densate and Euler-Heisenberg type effective action with 
quarks in loops. We develop a self-consistent scheme 
based on a full set of DSE with dynamical quarks, ghosts 
and gluons. Our approach is built on methods and pre- 
scriptions we adapt from Fischer 

Non-local quark condensate was first introduced by 
Mikhailov and Radyushkin 0] , and further developments 
followed soon after 0, 0]. The gauge invariant NLC is 
defined by 



C(x 2 ) = (q(x)E(x;0)q(0)), 
where Wilson phase factor is defined as 



E(x;0) = Pexp ( ie J A^(x)dx^ 



(1) 



(2) 



and the contour C connects points x and 0. In this pa- 
per we focus our attention on the first terms in powers 
of x which are independent of the Wilson line contri- 
butions. Wilson line terms are in general very impor- 
tant, and such contributions should be evaluated self- 
consistently @, which may be done by the Ericson- 
Semenoff-Szabo-Zarembo (ESSZ) technique used by us 
in the previous work It has been shown within the 
instanton vacuum model |7j that the form of the NLC is 
nearly independent on possible irregularities of the path, 
such as a cusp and thus in general the path can be rep- 
resented by a straight line. 

The initial motivation for introducing a NLC came 
from its influence on hadron phenomenology. For this 
reason NLC has been decomposed into the local conden- 
sates (LC) and the measure of the quark fluctuations in 



vacuum, known as the quark virtuality (QV). This quan- 
tity related to NLC, is defined as 



, _ (qDlq) 



K = 



(3) 



(here is the covariant derivative) , arising in the stan- 
dard operator product expansion (OPE) of the NLC as 
the coefficient in front of the quadratic term: 



C(x 2 ) = (g(O)g(O)) 



1 



(qD 2 q) 
(qq) 



(4) 



Quark virtuality is related to the gluon-quark trilinear 
(local) condensate 



(qD 2 q) 



W) 



(qga^G^q). 



(5) 



and thus can be counted as an independent vacuum 
structure parameter, characterizing the non-perturbative 
QCD vacuum. The standard estimate for X 2 by Chernyak 
and Zhitnitsky [8] is X 2 . » 0.4 ± O.lGeV 2 . There are 
larger estimates however, e.g. Shuryak suggests [9] 
X 2 ~ 1.2GeV . We note that these numerical values for 
the correlation length are comparable with the typical 
hadronic scale. 

Our effort to relate DSE and NLC is not the first. An 
attempt to derive self-consistent equations upon conden- 
sates was made by Pauchy Hwang [T3| in the large- 
l/N c limit; unfortunately, this was not developed fur- 
ther. Non-local quark condensate has been obtained 
within the flat-bottom potential approach to Dyson- 
Schwinger equations, where typical correlation length 
of SGeV" 1 has been obtained [ll[. Dyson-Schwinger 
equations are solved in [13] for quark dynamical mass 
and wave-function (no gluons or ghosts solved dynam- 
ically; gluon propagator mimicked by an Ansatz, rain- 
bow approximation applied to quark equations); using 
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propagators, the quark-quark non-local condensate and 
the quark-quark-gluon local condensates are calculated, 
typical correlation length obtained is O.SGeV" 1 . Same 
methods were used in 1131. where virtualities A 2 = 



0.7GeV 2 ,A 2 = 
confirmed in 



1.6GeV are reported. These results are 
[T3 | and completed with gluon virtuality 



as well A 2 = 0.2a7 5 - l.Oc^GeV , the latter exhibiting 
strong scale dependence via coupling constant. Dyson- 
Schwinger equations were solved in a similar approxi- 
mation (no dynamical gluons and ghosts) in [|l5j; how- 
ever, surprisingly large values of virtualities have been 
reported: A 2 d = 12 . . . 16GeV 2 , A 2 = 14 . . . 18GeV 2 . Till 
now, there has been no self-consistent treatment of the 
non-local condensates based on DSE with gluons. We 
consider this to be a disadvantage of the scheme, since 
quark fluctuations in vacuum are driven by gluons. We 
will present our result for QV as function of quark mass. 

We describe DSE methodology and calculate the prop- 
agators in the next section (JTTJ) , the non-local condensate 
(NLC) and its response to an external field is studied in 
section llIIII) . In section l|IVp we do the Euler-Heisenberg 
type effective action for quarks with non-perturbative 
DSE propagators in external fields and compare our re- 
sults to the meson based evaluation. We conclude in 
section (TV]) . 



II. DYSON-SCHWINGER EQUATIONS 
A. Formulation of DSE with Quarks 

In this section we review the technique of obtaining 
quark and gluon propagators in a self-consistent way. We 
use Fisher's DSE technique described in [3], and we show 
that the propagators are reproduced by us in the case 
with quarks as well - we have already reproduced the 
gluodynamics sector in our previous study 0]. 

We apply in our work the Newton optimization 
method, based on the numerical procedure described 
in [lB|. We solve a system for ghost, gluon and quark 
propagators, as shown in Fig. {TJ. Propagator dressing 
is shown by bulbs, and that of vertices - by transparent 
bulbs. We parameterize the gluon propagator in Landau 
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FIG. 1: Diagrammatic representation of DSE. 



gauge by the form-factor F, defined via relation 

P^pA F(p 2 ) 



•\F ab { ^\ ^ab 



D* ab (p) 



9^ 



P 



p z + it 



(6) 



and the ghost propagator is parameterized by the form- 
factor G 



D 



Gab 



(P) 



§ ab 



p" + ie 

Quark propagator is defined as 

' S {P) = 1 1 



-G(j? 



(7) 



(8) 



A(p) j> + M(p) ' 

Finding the scalar form-factors F, G, A, M will yield non 
perturbative information on the physical quarks and glu 
ons. 

DSE for this system can be written in the form: 



1 



1 



F{p 2 ) 
1 

Mx) 



A(x) 



(iLy)-n(^)), 
MW/i 2 ) 



(9) 



[ M(x)A(x) = M(fl 2 g ) + Jl M (x) - n M (f g ) 

Here fig iC the are points of subtraction, /2 C = 0, ft g = p,, ft 
is the limit of the interval p 2 € (0, ft 2 ) in the momentum 
space where we solve DSE, coupling g 2 is meant to be 
taken at point /i: g 2 {ft 2 ). Gluon vacuum polarization is 

n( P 2 ) = n 2 V) + n 2 V) + n 2 V), (10) 

contributions of ghosts II 2c (p 2 ) and gluons Ii 2g {p 2 ) being 
d d q 



Tl 2 %p 2 )=N c g 2 



(27T) 

Tl 2 °(p 2 )=N c g 2 J J0- iQo{ p 2 , q 2 ,r 2 ) F { q 2 )F(r 



M (p 2 ,q 2 y)G{q 2 )G(r 2 l 



^host self-energy is 

E(p 2 ) = N c g 2 



K (p 2 , q 2 y)G( q 2 )F(r 2 )^L. 



(11) 
(12) 



Quark self-energy is conveniently split into functions Ha 
and Hm> given below: 

1 f 4 f a(z) G{z)- 2d - d / s 1 
M "3^7 y \z(y + M 2 (y)) F(z) d A{y) 



~ (A(x) + A(y)) M(y) + \ (AA(x, y)M{y)- 
-AB(x, y)) (-z + 2(x + y)-(x- y) 2 / z) + 

&A(x)-A(y))M{v)Cl(x,v)(x-y) 



(13) 
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and 
Ha 



1 

3^3 



d 4 y 



a(z) G(z) 



-2d-o 



-z + 



x + y 



xz(y + M 2 {y)) F(z) d A{y) 

(x-y) 2 \ A(x)+A(y) 
2z 



^A(x,y) 



(x + y) + AB(x,y)M(y) x 



(_! + (» + „) + 



(A(x)-A(y))n(x,y) 



2 2 

x —y x — y 
~ z — - — 



(14) 

yielding the last two equations of ((9|). Here auxiliary 
functions AA, AB, Q,, Ail have been introduced: 



AA(x,y) 

B(x) 
AB(x,y) 



A(x)-A(y) 
x-y 

M(x)A(x), 

B(x) - B(y) 
x-y 



(15) 



(x - y) 2 + {M 2 {x) + M 2 {y)f 



The constructions lfT3|) . lfT4|) . are taken from [l(, we have 
fixed here a typo originally present in Eq. (fl3"ll . Parame- 
ter d is related to the Ansatz for the quark-gluon vertex 
that is used. There is no unambiguous way of choos- 
ing this parameter, since there is no fully consistent way 
of truncating DSE without violating some of the worthy 
properties of the original full tower of equations, and we 
refer the reader to [l| for a comprehensive discussion on 
that point. Variable z is a logarithmic variable 



z = In ■ 



(16) 



and scale fi is yet to be defined upon solving DSE from 
comparing the obtained coupling ocdse{z) to the known 
values of Particle Data Group coupling apDcip 2 ) [13] at 
point M: 



aosEiHM 2 /^)) = a PDG (M 2 ). 



(17) 



The coupling constant <? 2 /47r = a is expressed in terms 
of G, F solely [ill Ell, as vertex is finite in Landau gauge 
(at one-loop level) 



<*DSE(ln( P 2 )) = a DSE (/2)F(p 2 )G 2 (p 2 ). 



(18) 



In our case we shall use varying scale fixing point M so 
that we can prove that our results are independent of 
scale fixing point choice within the error margin of our 
procedure. 



The kernels Mo,Kq,Qo are well-known in literature, 
however to make the presentation self-contained we pro- 
vide them here: 



Ko(x,y,0) = 
M (x,y,9) = 

Qo(x,y,0) = 



y 2 sin 4 (9) 



(—2cos(9)^/xy + x + y) 

y 2 sin 4 (6>) 

3a; (—2 cos(9)y/xy + x + y) ' 



(19) 



12a; (-2 cos(9)^/xy + x + y) 
{y sin 2 (6>) [2cos(26>) (6a; 2 + 31a;y + 6y 2 ) - 

-12a; cos(39)y/xy + xycos(49) - A&co$(6)-^xy(x + y)- 

-12ycos(39)^xy + 3x 2 + 27xy + 3y 2 ]} . 

(20) 

Scalar variables x = p 2 , y = q 2 are introduced; variable 9 
is defined via (p—q) 2 = x+y — 2^xycos9. We neglect the 
effects of non-trivial dressing of the vertices, since these 
do not essentially back-react the upon the IR structure 
of the propagators themselves. 

To solve the Dyson-Schwinger equations we use the 
Ansatz (HEDl: " 



F(z) 



exp ^aiTi(z) \,z€ (lnejn, 



2 \ 7 

F(fi) I l + wlog^- ] ,z>ln/x 2 , 



Az 2k , z < lne 



(21) 



G(z) 



exp I y^bjTi(z) J , z E (lne,ln^ 2 



G(a) 1 + wln- 



, z > In /j 2 



Bz K , z < e. 



and similar Ansatze for M{p) 1 A(p). Here Tj are 
Tschebyschev polynomials, a,, bi are unknown coefficients 
yet to be determined from the numerical solution, n is the 
number of polynomials used (mostly n = 30 has been 
used here, allowing precision of up to 10 -10 for the coef- 
ficients), S = -9/44, 7 = -1 - 26, cj = llN c a(ay(Un). 
The IR scaling k is chosen to be the standard [22, 23] 



k = 0.59 



(22) 



for the case of Brown-Pennington truncation with £ = 
1 [2l| (for discussion of meaning of £ see [13]), which 
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is our case (C already set to its number value every- 
where). Following [25j |. we employ renormalization con- 
stant Z\ redefinition, so that no momentum dependence 
could possibly enter it, that is 



1 I A 



G(y) 



(1— a/6-2a) 



G(y) 



(\-b/S-2b) 



Again, following [H[ we choose 

a = b = 3(5, 



(23) 



(24) 



which minimizes its momentum dependence. Renormal- 
ization constant Z\ refers to the piece with a ghost loop 
in vacuum polarization. The equations are solved by us- 
ing Newton's method, described for this particular ap- 
plication by Bloch [16| . The results of the solution are 
propagator form factors F, G, shown in Fig. (J2j) on the 
left, the IR behavior of the propagators corresponds to 
the standard ghost enhancement and gluon suppression. 
The coupling a obtained from DSE (fT8|) is shown on the 
right in Fig. We note here that the IR fixed point 
seen in the Figure is 



«(0) 



(25) 



for iV c = 3, which is consistent with the up-to-date 
Dyson-Schwinger results reported by other groups (2H . 




0.1 0.2 



FIG. 2: Ghost (dashed line) and gluon (solid line) 
propagator form factors obtained in DSE in Landau 
gauge; running coupling from DSE. 

Quark wave-functions were obtained for one quark at 
a time solving a in a selfconsistent way the DSE, i.e. 
these are unquenched quarks. They are quite similar 
to quenched approximation where quark DSE is solved 
for given glue DSE solution (quenched approximation) 
The wave function form factors are shown in Fig. (fS]). 
Wave-function form factors become perturbatively unity; 
within an error margin they are no more distinguishable 
in the UV, although they exhibit different and non-trivial 
behavior in the IR. 

The quark masses are shown in Fig. ((4]). Physically 
it is important that UV anomalous dimensions of all the 
quarks are rendered the same in Fig. (@1) , which con- 
firms validity of the procedure. This can be seen from the 
dashed parallel lines in Fig. Q In general, in this Sec- 
tion, we confirm all the current knowledge on the DSE 




— /V 

— d 

— s 



p , GeV 



FIG. 3: Wave-function (propagator) form factors A(p) 
for flavors u(red), (i(green), s(blue), c(magenta) (lines 
from top to bottom at p — 0.5 GeV. 



with quarks. We improve the numerical convergence by 
smoothing numerical cut-off on integrals by superposing 
varying limits, which procedure removes Fourier trans- 
form 'echos' from the results. 

M, GeV 
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p . GeV 

0.15 0.20 0.30 0.50 0.70 1.00 1.50 2.00 3.00 

FIG. 4: Quark mass M(p) for flavors u, d, s, c. 
Punctured parallel tangent lines demonstrate that 
anomalous dimension is mass-independent. 



III. NON-LOCAL CONDENSATE 
A. Dependence on mass of condensate shape 

In this Section we calculate the non-local condensate 
omitting the Wilson line, study its behavior under ex- 
ternal fields and compute the vacuum response due to 
presence of non-local condensates to external fields. 

The nonlocal condensate-related vacuum expectation 
value (C-VEV) 

C (x) = $(x)1>(0)), (26) 

where local condensate C(0) satisfies C(0) = C (0), can 
be related to propagator as 



C (x) 



(20 



i Nf r 



e i P x 4M,(p) 
Mp) pi + M? 



(PT) (27) 
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However, separation the perturbative (PT) part from the 
non-perturbative propagator is not well-defined. More- 
over, some argue that the non-perturbative procedure is 
producing only the non-perturbative quark propagator 
and there is no PT subtraction needed. We do not have 
a good argument to support this reasoning, or, alterna- 
tively, a PT part subtraction, thus we follow the former 
approach. This also does not introduce additional pro- 
cedure ambiguity. Accordingly, it should be remembered 
when evaluating our results that the full non-perturbative 
understanding of QCD vacuum cannot be reached on 
grounds of Dyson-Schwinger equations alone, without ap- 
plying additional resummation procedures, e.g. ESSZ- 
resummation 0|. For this reason our results should be 
treated as a first qualitative estimate, and not yet as ex- 
act predictions. 

We think that despite any of the above shortcomings 
the results we obtain are surprizing. We show the C- 
VEV in Fig. |(5]), where we see from top to bottom 
(at x —> 0) beginning with the heavy quark (c(x)c(O)), 
(s(x)s(O)) (d(x)d(O)), (u(x)u(Q)), and last the lightest u 
quark. Numerical difficulties prevent us from reaching 
higher mass than 500 MeV for charm (at scale of 2 GeV) . 
Non-local condensate exhibits some oscillatory behaviour 
within 2 < x < 10 GeV" 1 . We believe that these C-VEV 
oscillations are due to numeric uncertainty, but their per- 
sistence and appearance when scale matching of quarks 
to glue and ghost occurs suggest that further study of 
this phenomenon is needed, and thus we show these re- 
sults with C-VEV reaching to 2fm distance. At this large 
distance the sequence of the C-VEV has reversed with 
smallest quark mass leading to largest value of C-VEV. 




FIG. 5: Non-local condensate for u, d, c, s quarks. 



B. Local quark condensate and quark virtuality 
dependence on mass 

The standard wisdom [13] about condensate depen- 
dence on mass for heavy quarks is 



(qq) = 



1 



This relation is usually derived from requiring continuity 
between heavy and light quarks' properties, imposed at 
the scale of about 0.2 GeV. The behavior of our propaga- 
tors and wave functions is continuous, yet the dependence 
on mass we observe is completely different. Another well 
regarded relation is [28] 



iss 



0.8(uu) 



(29) 



Note that in our evaluation the local condensate is in- 
dependent of Wilson line integral and thus our results 
for x — > while still PT subtraction dependent are more 
secure. For c and s quarks the values one sees in Fig. |(5]). 
are considerably larger than expected. Moreover, we find 
that our condensates increase with mass and does not de- 
crease, as was expected based on above estimates. 

The values of condensate is fitted surprisingly well by 
a simple power law 



n / rn \ 0.73 

o(m) = 0.2GeV»(_) 



(30) 



not at all expected from any qualitative QCD model 
we know. Mass dependence of condensate is illustrated 
in Fig. ([6]). The dashed line is the expected light quark 
value, the thick line the c, s expectations of Eq. (|28|) . Note 
that these results are obtained by considering one quark 
at a time and solving selfconsistently DSE (unquenched 
single quarks). 
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FIG. 6: Local condensate mass dependence, red dots 
are DSE results, thick line - estimate lf2"8| for heavy 
quarks, dashed line - standard (uu) value, thin line - 
power-law approximation (|30| . 

It seems that with increasing mass quarks can probe 
better the non-local glue vacuum fluctuations and thus 
their response strength increases. The non- locality of 
glue vacuum structure is usually not considered in the 
qualitative condensate models. However, there is no ar- 
gument we can present to align the light quark local con- 
densate as function of m with heavy quark condensate. 
We also note that when dealing with realistic quarks, 
their physical magnetic moments must be taken into ac- 
count. However this effect diminishes with quark mass 
and cannot explain the heavy quark condensate behavior. 
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In another attempt to understand this strange behav- 
ior, one could suggest that heavy quarks are worse repre- 
sented by Dyson-Schwinger equations since they tend to 
decouple and thus one-loop approximation becomes al- 
most free, but at higher loops they might become again 
important, thus yielding DSE approach invalid. How- 
ever, this explanation is not valid, since comparison of 
quenched approximation to the unquenched shows very 
little difference between the two. Thus the issue of con- 
densate dependence on mass in DSE scheme presented 
here remains an open question. 

Should this behavior be true, this strong dependence 
on mass of light quark condensate would deeply impact 
the chiral model analysis of quark masses, where a cor- 
nerstone assumption is that light quark condensates have 
equal value. 

Quark virtuality dependence on mass is given in the 
table U below and is shown in Fig. ((7]). We note the 
highly regular behavior, following the fit 



local condensate 



A 2 = 0.39GeV 



( m i \ 
VlGcV/ 



1.07 



(31) 



shows in Fig. J7]) For comparison recall that virtualities 
X 2 ud = 0.7GeV 2 ,\ 2 s = 1.6GeV 2 were reported [3], as 



discussed in Section 1. Recall also A 2 
and A 2 ~ 1.2GeV 2 [|. 



X 2 n , GeV 2 



0.4±0.1GeV^ 




0.002 
0.001 



FIG. 7: Quark virtuality dependence on quark mass. 



d 2 



W)f = W)o + F 2 h{x) - F ua F« — f 2 (x), (32) 



where moments /i, fa are 



h(x) 



1 r e ipx d 4 p (-8)m(p) 
(X) 1 J A 3 {p) (p 2 +m 2 (p)) 3 ' 

1 r e ipx d 4 p (-16)m(p) 
(2^0* 7 A3{p) {p 2 + m 2 {p)y 



(33) 



We notice here that not only the character of condensate 
dependence on x changes due to field switch-on, but it 
acquires anisotropy. The function /i is shown in Fig. |(9|). 
It deserves attention that smallest quark masses bring 
largest response to field, which is quite reasonable. The 
resulting parameters are shown in the table |l]). It can 
be seen from analysis of /i that already fields of order of 
magnitude of 10 -1 GeV 2 may put the local condensate 
to zero. This is comparable to the prediction of critical 
fields for condensate 



log 2 



(34) 



by Smilga and Shushpanov [2£ 



Si//i///S(eF ) , GeV 




x , GeV 



FIG. 9: Factor fi(x) as function of distance, describing 
the nonlocality of condensate sensitivity to field. 



C. Condensate response to external field 

We can also establish the character of condensate 
dependence on the external field. Considering a dia- 
gram Fig. ©, we derive the F 2 -order term in the non- 




FIG. 8: Diagram describing condensate sensitivity to 
field. 



Long-distance correlations will be even more sensitive 
to fields, since /i decreases slower than condensate it- 
self, thus making pion wave-function a nice candidate for 
analysis in an external field. 



IV. EFFECTIVE ACTION DUE TO 
CONDENSATES 

One of the simplest nonlinear processes of QED is 
photon-photon scattering, shown in Fig. (|10)| . In the lan- 
guage of Euler-Heisenberg effective action, the following 
term is responsible for this kind of processes 



C = a{F^F^) 2 + F V ,F^F P F^ 



li v X p 



(35) 



A{F„ V F^) 2 + B(F^F^) 2 
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TABLE I: Main characteristics of condensate: local amplitude C(0), virtuality A , condensate amplitude variation 

S ?pi) , virtuality variation jj$z, infrared exponent a ((^p(x)^p(O)) ~ e~ ax ), variation of the infrared exponent jp?. 
Mass value m = 0.51 in the fourth line is not a misprint against the expected m — 1.27 GeV, but was the largest 

mass at 2 GeV scale available to us. 



m q , GeV 


q 


C(0),GeV 


A 2 , GeV 2 


^GeV" 1 


% GeV" 2 


a, GeV 


-^r GeV -3 


0.0025 


2/3 


0.00239 


0.00066 


0.037 


0.00082 


0.40 


0.20 


0.005 


1/3 


0.0042 


0.0013 


0.023 


0.00094 


0.65 


0.20 


0.105 


1/3 


0.037 


0.039 


0.015 


0.0017 


1.04 


0.22 


0.51 


2/3 


0.12 


0.18 


0.0064 


0.0015 


1.04 


0.37 



Coefficients a, b are in case of QED 



b = 



36m 4 ' 
la 2 



(36) 



90m 4 



and A, B are linearly related to them: A = a+b/2, B = b. 
We shall calculate now these coefficients for condensate 
contribution of QCD vacuum into QCD-related photon- 
photon scattering. We shall see at the end that the con- 
tribution is larger than expected, compared to standard 
(perturbative) contribution due to hadrons. However, 
the magnitude of the effects we find is very small com- 
pared to what is experimentally accessible today, and 
in the foreseeable future, in the domain of intense laser 
physics. 



41 / 
+k, 



K jq 
\ 



p+K 



p+k,+k 2 



FIG. 10: Leading nonlinear term in the 
Euler-Heisenberg effective action, ki are incoming 
momenta, q quark charges. 

Strictly speaking, when dealing with realistic quarks, 
their physical magnetic moments must be taken into 
account. In the effective action quark magnetic 
moments would invoke a contribution of the type 

^J?P-\V\ pllzV'j, plliVi ^ r J (T ^il/i (T /J2l'2 cr M3^3 (T M4l'4j > Nq^. 

ing that the outcome might be non-negligible, we leave 
this contribution aside, since it requires a serious modifi- 
cation of the DSE solution scheme and complicated issues 
of truncation validity. 



To achieve the result we calculate the diagram Fig. lflO|) 
with propagators obtained in the previous section, which 
are responsible for condensates. We do not separate 
the condensate and the free terms at the level of each 
propagator, but rather do the full diagram with the full 
propagators, and then compare to the perturbative terms 
(a, b already given above, multiplied by respective quark 
charges) . As momentum dependence of the full diagram 
on s, i, u invariants would be known only numerically as 
a result of a calculation, containing numerical data for 
propagators, we use the following trick. We work out the 
scattering amplitudes M(ei, e2, e3, cj), given as 



M^x P (kuk 2 ,k 3 , ki) = e 4 / ^ tr [S{p)^S{p + k{yf 
Sip+h + ^Sip + k^j"] . 

(37) 

For the scattering amplitudes at small values of photon 
frequencies u> we extract the coefficient at the uj 4 term: 



M^xpih, k 2 , k 3 , k^e^e^el = Af +w 4 a 2 M(ei, e 2 , e 3 , e 4 ), 

(38) 

for two specific sets of polarization vectors, namely, 
(ei_L, e2_i_, e 3 j_, e 4 j_) and (ex||, e2\\, e 3 ±, e4±), (e;± denotes 
polarization orthogonal to reaction plane, and e x | polar- 
ization in reaction plane) , at specific values for 9 (namely, 
forward scattering 9 = ir) . These can be expressed as fol- 
lowing scalar integrals 



M(ei_L,e 2 x,e3x,e4j_) = 
32p 3 dp 



[I9p 8 + 75M(p) V 



o 15 [p 2 + M{p) 2 f A{pf 
10M(p) V - 330M(p) V + 30M(p) 



(39) 
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TABLE II: Coefficients a, b, A, B of non-linear terms in 
the effective action. The "perturbative" line shows for 
comparison the coefficients oo, fro, Ah Bo for mass m in 
the loops which can be thought of approximately as 
A ~ 300GeV; our results are shown as dimensionless 
ratios against ao,&o, Aq,Bq, Quarks charges 
qi = 2/3, 1/3 are included into coefficients. 





Oo 


bo 


A 


Bo 


PT 


1 


7 


i 

90mA 


7 

90mA 


flavor 


a/a 


b/b 


A/Aq 


B/B 


u 


0.07732 


0.093f7 


0.f328 


0.093f7 


d 


0.00302 


0.00337 


0.00425 


0.00337 


s 


0.000f9 


0.00022 


0.0003 


0.00022 


c 


0.00064 


0.0007 


0.00085 


0.0007 



Basing on the analysis of light quark properties we 
confirm the claims that QCD vacuum cannot be probed 
via non-linear QED effects until fields of hadron scales 
to 2 ~ (O.lGeV) 2 are reached, actually on scales several 
times below typical hadronic scale the effect can be felt 
as contribution to the naive QED estimate of photon- 
photon scattering. 

We believe that the issue of heavy quarks, which have 
yielded several inconsistencies in the suggested scheme, 
could be addressed in our scheme once more effort is de- 
voted to understand the solutions of DSE. Moreover, to 
arrive at realistic results we need to consider schemes 
with all quarks participating in the solution of DSE equa- 
tion, and allowing for vertex correction which implements 
magnetic moment for light quarks. 



V. CONCLUSION 



M(e 1 ||,e 2 ||,e 3 x,e4_L) = 

Jo 15 [p 2 + M(p) 2 f A{pf 1 
-40M(p)V + 60M (p) V ~ 30A/ (p) 8 ] . 

(40) 

Polarization vectors have been 

e\\ ={0,0,1,0}, 
e ± = {0,0,0,l}, 

with center-of-mass kinematics 

fci =cj{1, 1,0,0}, 
ka=u{l, -1,0,0}, 

= ui{1, cos9, sin (9, 0}, 
A)4 — ll>{1, — cos 9, — sin 6, 0}. 

In the expansion we used the fact that to is believed 
to be small, therefore, all non-perturbative momentum- 
dependent factors (M(p),A(p)) are taken at the point 
P- 

On the other hand, the coefficients M(. . . ) are known 
from ([35)) by direct analysis 

Af(eix, e2_L, eax, ea) = 64(2a + b), 

(43) 

-W(ei||,e 2 ||,e3x,e4x) = 16(4a + b). 

Thus a simple comparison of ([39]) and lf43|) yields val- 
ues for a, b and A, B. They are shown in Table pi)) . 
This Table is quite instructive. First of all, the con- 
tributions are comparable with the expected hadronic 
ones. The range of the latter can be estimated roughly 
within l/90m* . . . l/90m^ - 40 . . .0.05 GeV" 4 . Quarks 
with large bare masses yield less, as expected on general 
grounds. 



Solving one quark-gluon-ghost Dyson-Schwinger equa- 
tions, we have obtained the quark non-local condensate 
and quark virtuality as functions of quark mass. The 
mass dependence of the condensate disagrees with cur- 
rent qualitative wisdom. We could not find an explana- 
tion for why this is the case. The growth of the quark 
condensate with m 73 ,m < 500 MeV implies a signifi- 
cant difference between all mass condensates above and 
beyond any expectations. If confirmed, this result would 
have considerable impact on hadron phenomenology. For 
example, the difference in u-quark and (i-quark mass and 
thus implied difference in quark condensate leads to con- 
siderable change in chiral analysis of quark masses. 

Regarding the influence of an external field on the con- 
densate we predict that fields of order of magnitude of or- 
der of magnitude of 10 _1 GeV 2 can actually destroy local 
condensate, and even smaller fields can destroy non-local 
x-dependent condensate at x ^ 0. This result may have 
direct impact the pion wave function in external fields. 
As recently shown by Pimikov, Bakulev and Stefanis [30| . 
it is the non-locality of the condensates that is the key 
point for inclusion of the non-perturbative contributions 
to the pion form factor. Thus our results on condensates 
immediately drive the dynamics of pion wave function in 
external fields. 

In addition we predict that light quarks yield im- 
portant non-perturbative contributions into the photon- 
photon scattering amplitude, which are comparable with 
the corresponding perturbative contributions based on 
loops with light mesons. We stress that this dynam- 
ics is essentially condensate-driven, and particularly by 
its non-locality. Our present non-perturbative evaluation 
suggests that the critical field, above which the non-linear 
QCD-QED effects can be seen, is several times lower than 
the typical hadronic scale. Even so, the experiments to 
probe the QCD vacuum with intense laser fields are be- 
yond the foreseeable future. 

We have outlined in the text an opportunity for fur- 
ther advance which must first focus on the resolution of 
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the mass dependence of quark condensate and better un- 
derstanding of the related virtuality. The relatively large 
effects which external fields can impart on the QCD vac- 
uum must be confirmed in the context of such an im- 
proved theoretical framework. 
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